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I. INTRODUCTION 

In recent years, branes have been studied due to their 
uses in diverse topics of high energy physics. Those more 
relevant are: string dualities, which ended up with the 
breaking idea of the AdS/CFT correspondence [1, 2], and 
a new interpretation of compactification, which allows to 
postulate the ^-theory [3, 4] and the realization of brane- 
world scenarios [5, 6]. 

These branes are extended objects arising in super- 
string theories, but their low energy analog might by real- 
ized through a gravitational domain wall. Gravitational 
domain walls are solutions to the Einstein-Hilbert-Klein- 
Gordon equations where the self-interacting term of the 
scalar field has (at least) two degenerated minima and 
the scalar field interpolates between consecutive minima. 

It has been shown that gravity does get localized on a 
D3-brane embedded in a higher dimensional space-time, 
in the sense that the effective theory of gravity on the 
domain wall is described by a Newtonian potential plus 
corrections, which coincide with the first one due to gen- 
eral relativity. However, in order to provide a good phe- 
nomenological scenario, they should localize fermions and 
gauge bosons as well. Several generalizations to these 
models have been developed, such as considering more 
than two D3-branes[7], domain walls that are the thick 
versions of branes [8], domain walls generated by more 
than a single scalar field [9], and so on. 

If thin 3-branes are considered, fermions cannot be lo- 
calized [10] and therefore, generalized brane-world sce- 
narios have to be analyzed. Among the possibilities one 
might generalize the brane to their thick analog. 

When a domain wall is considered, the simplest new 
interaction to be added is a Yukawa-like interaction be- 
tween the scalar field which generates the domain wall 
and the fermions added to the formalism. However, in 
generic models only a single chirality is localized. In a 
previous work [11], the moduli space of parameter for a 
particular family of domain walls, known in the literature 
as double domain walls, was studied. It was found that 



there exist regions where possibly both chiralities are lo- 
calized, and that the most plausible scenario was highly 
dependent on the Yukawa-like coupling constant. 

The aim of this work is to provide a similar analysis 
of the moduli space for diverse kinds of domain walls, 
which is an essential ingredient in order to find out those 
that are acceptable candidates for phenomenology. In 
the next section a set of different domain wall solutions 
will be introduced. In section III a brief review about 
localization of fermions on static branes is considered. 
In the subsequent sections analysis and discussion are 
presented. Finally, the notation is include in appendix 
A. 



II. GRAVITATIONAL DOMAIN WALL 
SOLUTIONS 

Gravitational domain walls are solutions to the 
Einstein-Hilbert equations coupled with a self-interacting 
scalar field, 
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such that the potential of self-interaction has (at least) 
a couple of degenerated minima, and the scalar field in- 
terpolates between two consecutive minima. Since one 
is interested in solutions with a plane-parallel symmetry, 
the scalar field depends only on the extra-dimension, £ 
[12]. 

It is well known that the most general metric satisfying 
these conditions is, 
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(2) 

Static and dynamical configurations can be found by con- 
sidering B as B(£) or B(t) respectively. A bunch of gravi- 
tational domain wall solution might be found in the liter- 
ature, and they can be generalized to various dimensions 
[13]. Additionally, via scaling a thickness parameter can 
be introduced into the solutions, and families of solutions 
might be found [14]. 

In the following the restriction to five dimensions is 
assumed, and only static solutions are considered. 
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A. Kink wall 

The kink domain wall is the simplest generalization of 
the solitonic ^2-kink found in a 1+1-dimensional model 
[12, 15]. 

This solution can be achieved by setting in (2), 



A(£) = - -5 In ^cosh 
B(0 =C(0 = 0, 
and scalar field, 
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where 5 parametrizes the thickness of the wall, and a 
is related to the cosmological constant, which for this 
model is negative. It will be shown in the next section 
that the self-interacting potential of the scalar field is 
not important for the purpose of fermion localization, 
and therefore its explicit form is omitted. 



B. Sine-Gordon wall 

The sine- Gordon domain wall is also the generalization 
of the 1+1-dimensional sine-Gordon model [12]. Since 
the self-interacting potential has an infinite number of de- 
generated vacua, there are different static solutions which 
are absent in a Z2-model, such as the kink-kink solution. 

In this case, the functions in (2) are 



A(£) = -S\n ^cosh ^ 
B(0 =C(f) = 0, 
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This is the single kink-like solution, where the scalar field 
interpolates between a couple of consecutive minima of 
the potential. 



Nevertheless, an example of asymmetric domain wall 
which does localize gravity has been found [17], with 

A(0 = <*£ - S exp (-2e~ # ) + S Ei (-2e~ # ) , (9) 
B(0 = C(0 = o, 

and scalar field, 

0(0 = 2V35 (exp (-e - ^) - ej . (10) 

Here the function Ei(£) is the exponential integral, de- 
fined by 
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and e is a shifting constant, which determines the asymp- 
totic behavior of the fermionic modes. 



D. Double walls 

Double domain walls are a parametric family of so- 
lutions, whose energy density peaks in two £ =constant 
regions inside the wall thickness, and are a slight general- 
ization of the thick brane of the Randall- Sundrum model 
presented in [8]. These solutions are given by 
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and 



B(0 = 0, 



together with, 
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In order for cf) to interpolate between to minima of the 
potential, s must be odd. 

Fermion localization on this family of solutions has 
been studied in a previous paper. Readers interested 
in a more exhaustive analysis of these domain walls are 
referred to [11]. 



C. Asymmetric wall 

Although it has been shown that most asymmetric do- 
main walls do not localize gravity [16, 17], it is known 
that they could emulate a gravitational attraction be- 
tween particles on the wall. Thus, localization of fermions 
(and gauge fields lately) could result in a model where 
the fundamental content of the standard model lies on 
the brane, whilst four dimensional gravity is an effect of 
the embedding geometry of the brane [18]. 



III. FERMION LOCALIZATION ON STATIC 
D-BRANES 

In the following, thick-branes (or gravitational domain 
walls) are considered. The formal relation between them 
might be established through a distributional limit [14]. 

Using the equations of motion coming from the five 
dimensional Dirac Lagrangian on curved space-time with 
a Yukawa interaction, 

(- y-M + AP(0))^ = O, (13) 
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and the decomposition of the five dimensional spinor in 
terms of four dimensional chiral ones, one gets the profile 
equations [19, 20]. These equations are Schrodinger-like, 
and their quantum mechanical potentials depend on the 
four dimensional chiral condition and Dirac equation. In 
the following, a brief review of the conditions that should 
be satisfied for (chiral) fermions to localize are shown, but 
a detailed analysis is presented in [11]. We are only going 
to consider flat domain walls, i.e., B(£ = 0) = 0. 

A. Massless 5D Fermions 

When four dimensional chiral fermions, -0±, satisfy, 

^ (4 V± = o, 

the profile equations are (see for example [19, 20]), 

f± + A'f± T W{4>)e c f ± = 0, (14) 
where f± are the profiles for ijj±. Their solutions are, 

/ ± oce-^ ±A / d «^ eC . (15) 

However, four dimensional massive fermions mix chi- 
ralities, and the profile equations are coupled. Through 
a series of changes of variables, £' = J d^e~ A+c and 
f± i — y e~ 2A u±„ the equations might be decoupled, re- 
sulting in a Schodinger-like equation, 

[-dl,+V q ± m ]u ± =m 2 u ± , (16) 

where 

f±^e- 2A u±, 

and 

= (WWe A ) 2 ± % {W{4>)e A ) . (17) 



IV. ANALYSIS 

One might constraint the parameter space using ana- 
lytic methods. In the following, the polynomial function, 
V(4>), will be restricted to a monomial, 

V(<P) = 4> n , (20) 

to avoid the introduction of more free parameters. 

In order to assure that one of the chiralities is localized, 
a constraint is found for the coupling constant, say, 

/±(0 oc exp (-A ± J d^\V^)e c ^ , (21) 

when M = 0. Additionally, in the simplest model a small 
coupling constant value does not allow localization of the 
positive chirality, independently of the domain wall con- 
sidered. Thus, if both chiralities are expected to be lo- 
calized, a new constraint on A is found. 

For the numerical analysis of the moduli space, a Nu- 
merov algorithm together with a bisection method is 
used for finding the eigenvalue of the ground state of 
the Schrodinger-like equation for the profiles. The pa- 
rameters are restricted to A G [1.00,5.00] sliced in steps 
of 0.05, and M e [0.00,1.00] in steps of 0.02. These 
algorithms were implemented using SAGE [24]. 

The existence of tachyonic modes has been used to rule 
out several models, since they are a signal of instabilities 
in the model, due to causality. Nonetheless the recent 
findings of OPERA [25] could be interpreted as possible 
evidence of tachyonic behavior of neutrinos through mat- 
ter. In the figures, the vertical axis takes one of the two 
possible values, if the choice of parameters implies ex- 
istence of tachyonic instabilities, while 1 means that the 
choice of parameters is phenomenologically plausible, i.e. 
tachyon-free. 



B. Massive 5D Fermions 



Kink domain wall 



When a five dimensional fundamental Dirac mass, M, 
is added, it can be shown that the above equations, 
change by substituting 



XP((t)) XP((t)) - M. 



(18) 



From (15), it is clear that after the substitution (18), 
the localization of the positive profile, /+, is favored in 
detriment of the negative one, /_. This opposite effect 
serves as argument to suppose that regions (or a region) 
of the moduli space where both chiralities are localized 
exist (s). In the case of a massive fermion, the effect enters 
through the quantum mechanical potential (17), which is 
less intuitive to analyze. 

Nonetheless, it is well known that adding a fundamen- 
tal mass, M, in odd dimensional space-time gives a parity 
anomaly [21, 22], which vanishes when a sign function on 
(j> multiplies M[23], i.e. the substitution is, 



AP(0)^AP(0)-Msgn(0). 



(19) 



Using the warp factor and the scalar field correspond- 
ing to the kink domain wall, i.e. (3) and (5), in (21), one 
gets a lower bound constraint for the coupling constant, 
due to the negative chirality, 



A> 



2a 1 

Y (3S) n / 2 ' 



(22) 



This is the generalization of the constraints found in 
[20]. A similar constraint was found in [11] for the double 
wall scenario. 

In figure 1 it can be seen that the left-handed profile 
does localize itself, without the existence of tachyonic 
modes, for any value of M if A > A c ~ 1.44. On the 
other hand, the right-handed profiles are not localized in 
a tachyon free theory. 

Unlike the n — \ case, for cubic and fifth scalar- fermion 
coupling the left chirality does not present tachyonic 
modes while right profiles inevitable present tachyons, 
for any value of (A, M). 
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FIG. 1. Moduli space for kink domain walls for linear FIG. 2. Moduli space for asymmetric domain walls for linear 
Yukawa coupling and sign function, for left- and right-handed Yukawa coupling, for left- and right-handed fermions. 
fermions. 



B. Sine- Gordon domain wall 

Using (6) and (8) in (21), yields the constraint for the 
profile. However, integrating the argument of the expo- 
nential becomes highly non-trivial for monomials other 
than linear. Therefore, an approximation might be done, 



(23) 



Equation (23) defines the function T. Thus, the asymp- 
totic behavior of T determines whether or not the expo- 
nential converges. Convergence demands, 



A > a 
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The numerical analysis for this kind of domain wall 
assures that there is not a single point in the considered 
region of the moduli space which provides a tachyon free 
model of fermions. 



C. Asymmetric domain wall 

It was shown in [20] that for linear fermion-scalar cou- 
pling, (21) can be integrated explicitly. However, for 
more general couplings this is not the case, and (23) 
should be used. 

Due to the asymmetry, the limit values of A{£) and 
(/>(£) as £ — > ±oo differ, and conditions for localization 
vary for each limit and chirality. From (23) one get, 

-a + p ± A (2^35(1 - e)Y <0, (25) 
-a±X (2\/3fe) n <0, (26) 

for £ — >• +oo and £ — >• — oo respectively, which in terms 
of the cosmological constants and for negative chirality 





FIG. 3. Moduli space for asymmetric domain walls for cubic 
Yukawa coupling, for left- and right-handed fermions. 



are, 



A > 
A < 
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If one demands the saturation of both equations, a con- 
straint on £ is found, 



1 



A+ 
A- 



(29) 



and the critical value of the coupling turns out to be, 



1 



((A + )£ + (A-)*)' 



(30) 



The numerical analysis of the moduli space for this 
asymmetric domain walls gives a very interesting fea- 
ture. As can be seen in figures 2, 3 and 4, the left handed 
mode is either localized (for linear and cubic coupling) or 
unlocalized (fifth coupling). However, the right handed 
fermion presents a sort of "band" structure [26] , depend- 
ing on whether the model is tachyon-free (upper edge of 
the cube) or contains fermionic tachyonic modes. 

So far, the moduli space of these models has been sep- 
arated into pieces where tachyons can or cannot exist. 
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FIG. 4. Moduli space for asymmetric domain walls for fifth 
Yukawa coupling, for left- and right-handed fermions. 




FIG. 5. Moduli space for asymmetric domain walls for linear 
Yukawa coupling and sign function, for left- and right-handed 
fermions. 



When the sign function, which avoids a parity anomaly, 
is turned on, the left handed fermion presents always a 
tachyonic mode. On the other hand, the localization of 
the right handed profile for n = 1 has a band structure, 
where tachyonic modes can or cannot be avoided in the 
phenomenological model. For any other value of n no 
tachyon free model can be constructed. 



V. DISCUSSION 

Nowadays, due to the results of the OPERA experi- 
ment, physicists should ask themselves whether or not to 
model the universe allowing tachyonic modes. This new 
possibility opens the Pandora's box of causality, and of 
course new theories are been developed in this direction. 

Here, by considering brane- world scenarios with an ex- 
tra dimension, the simplest moduli space of parameters 
has been studied, and the numerical analysis has shown 
that the moduli space of parameters can be split into 
areas whose physical models are tachyon free and areas 
where this feature cannot be achieved. 

Several kinds of domain wall models were considered, 
and it seems that the most interesting structure of the 
moduli space is the one of the asymmetric domain wall, 
which presents what looks like a band structure, simi- 
lar to the electronic structure of condensed matter. Al- 



though higher fermion-scalar couplings were studied, the 
numerical results show that the most plausible scenario 
is the linear coupling among these fields. 

The sine- Gordon domain wall presents a moduli space 
which does not permit (in the considered region) to model 
physics without tachyonic fermions. Unless more data 
confirm the existence of the tachyonic modes in the neu- 
trino sector, the sine- Gordon domain wall scenarios must 
be ruled out for physical models. 

In kink domain walls, the n = 1,3,5 right handed 
fermions require the existence of tachyonic modes, while 
the left handed one is tachyon free, except for the case 
n = 1, which additionally requires A > A c ~ 1.44. There- 
fore, for A < A c , both fermion chiralities present tachyon 
modes. 

Some phenomenology can be read from the numerical 
analysis, despite the fact that the constraints on tachyon 
free models are strong. However, when models are con- 
structed it is important to have lines which allow to rule 
them out. 

Depending on the algebraic form of the quantum me- 
chanical potentials, V^, the localization of the left 
handed chirality is favored in detriment of the right 
handed one, due to the appearance of a barrier near 
the position of the D3-brane. Since light right handed 
fermions tend to propagate out of the brane, one might 
suggest that the OPERA results might be interpreted 
as propagation of right handed neutrinos through the 
warped extra dimension, as suggest in [27] as an expla- 
nation of the LSND neutrino anomaly. 
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Appendix A: Conventions and Notations 

Through the manuscript, the metrics have signature 
mostly positive. Since the tetrads formalism is used ex- 
tensively, distinction between fiat and curved coordinates 
by Latin and Greek indices. Moreover, hated indices run 
over the whole spacetime whilst unhated ones run over a 
hypersurface restriction, i.e., on the coordinates parallel 
to the topological defect. 

The gamma matrices are defined in the tangent space, 
and they satisfy the Clifford algebra, 

[r, 1 b }=2r 1 a h. (Al) 

In even dimensions one can define the chirality matrix 
7*, satisfying the properties 
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{ 7 a ,7*} =0. 



algebra are, 



(7*) 2 



From this, the projector operators, 



P± 



1-7* 



are both non-trivial. 



(A2) 



On the other hand, odd dimensional Clifford algebras 
are constructed by using the gamma matrices of the co- 
dimension one spacetime, via j a — (7 6 ,7*). These odd 
dimensional Clifford algebras have a trivial projectors 
P±, and therefore chiral fermions cannot be defined. 

In any dimension one may define a set of generators of 
the Lorentz algebra, constructed with the elements of the 
Clifford algebra (Al). These generators of the Lorentz 
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(A3) 



When the Dirac-Feynman slash notation is used it 
must be interpreted as, 



(A4) 



The covariant derivative for gauge theories is defined 



by, 



dp, - igAp 



(A5) 



with the gravitational connection, which is related to 
the Christoffel connection for integer spin fields, and with 
the spin connection for semi-integer spin fields. Clearly 
the Dirac-Feynman slash notation can be used with the 
covariant derivative, 

7= ^ 7 a V A = (0 A - igAf, - . (A6) 
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